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Abstract 

The conventional semi-empirical formula for collapse distance [Phys. Rev. 179, 862 (1969); 
Prog. Quant. Electr. 4, 35 (1975)] has been widely used in many applications. However, it is 
not applicable when the dispersion length is smaller than or has similar order-of-magnitude as the 
collapse distance. For the “enough short” pulses, there exists a threshold for the initial peak power, 
with which the collapse distance has a maximum value due to the competition between the Kerr 
self-focusing and the group velocity dispersion. New semi-empirical formulas are obtained for the 
collapse distance of the pulse with the initial power being less or larger than the threshold, and 
they can match the numerical simulations gracefully. 

PACS numbers: 42.65.Jx, 42.65.Tg 


* Corresponding author (email: wl@swjtu.edu.cn). 


1 




I. INTRODUCTION 


The ultra-short high intensity laser pulses have many applications such as laser-guided 
electric discharge [l|, terahertz generation ^ and remote sensing j^. Such applications need 
extremely high light powers which can be obtained via the collapse of pulses 

For a Gaussian beam, if the initial peak power exceeds the critical power for self-focusing, 
the beam will undergo collapse until the higher-order processes such as plasma or high-order 


n 


Kerr effects halt the collapse jSl-llOj . In most cases, the collapse distance Lc (the propagation 
length of the self-: 
formula 
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bcusing beam until collapse) can be well approximated by a semi-empirical 
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where uq is the refractive index, tq is the beam width which is at 1/e^ level of intensity, 

Ao is the laser wavelength in vacuum. Pin is the initial pulse’s peak power, and Per is the 

critical power for self-focusing which can be written as Per = 3.77AQ/87rnon2 with n 2 being 

the Kerr index. However, for the very high powers such as Pj„ = lOOPce, experiments and 

numerical simulations show that the collapse distance can not be described by and a 

transition from a 1 /\/Pj„ to a 1/Pin scaling was observed [Iq. Recently, simulations show 


that the group-velocity-dispersion (GVD) has a great influence on Lc when the pressure is 
relatively high, e.g., 10 atm jl^ . 

In this paper, we investigate the collapse distance of femtosecond laser pulses in air 
for different temporal durations. We find that for the “very short” pulses, there exists a 
threshold for the initial power, with which the collapse distance has a maximum value. More 
importantly, new semi-empirical formulas are obtained for the collapse distance of the pulse 
with the initial power being less than or larger than the threshold. The new formulas take 
into account of the influence of GVD and can match the numerical simulations perfectly, 
including the cases of high pressures at which the GVD becomes more important. 


II. THEORETICAL MODEL 

The physical processes which halt beam collapse is still controversial. Some groups think 
the plasma prevents collapse (Kerr-plasma model) 18l-l24|. while other groups believe the 
collapse is stopped by the high-order Kerr effects (HOKE model) 25l-l31|. Fortunately, 
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the collapse distances are almost the same for Kerr-plasma model and HOKE model 25, 
32|, since the HOKE and the plasma effect becomes important only after the pulse has 
collapsed 33j, and this phenomena has also been conhrmed in our numerical simulations 
(not shown here). In this work we study the collapse distance of the pulse for different 
temporal durations via HOKE model, which can be described by an extended non-linear 
Schrodinger equation (NLSE) 25, 27| 


dA _ i . . _ ik" d^A 

dz 2fco 2 (9r2 

i/Cp UJpe ^ A 


i=02,N2 


Wi{I)Ui 




'(^Pat,l Pe,l} 


( 2 ) 


where A represents the envelope of the electric held, and z denotes the propagation distance. 
ko = 2 t[/\q (Ap =800 nm) is the central wave number. The Laplacian operater = 
-|- 1/rdr denotes the beam transverse diffraction. The remaining terms in the right-hand- 
side of Eq.(l) account for the group velocity dispersion with the second order dispersion 
coefficient k” = 0.2 fs^/cm, Kerr and high-order Kerr effect with nonlinear refractive index 
77.2 = 1.2 X 10“^® cm^/W, ru = —1.5 x 10“^^ cm'^/W^, rip = 2.1 x lO”"^® cm®/W"^, and 


rig = —0.8 X 10 cm®/W^ 3J, [SSj, plasma defocusing with the plasma oscillation frequency 


ujpe = (ge is the electron charge, mg is the electron mass and p is the free electron 

density.), the energy loss caused by MPA. and Wo 2 {.I) are the photoionization rate 


of N 2 and O 2 . Pat,N 2 = 2.1 x 10 


25 j^-3 


and Pat^Oi = 5.7 x 10 


24 m 3 


are the density of N 2 


and O 2 molecular at 1 atm. Pe,V2 and Pe,02 are free electron density ionized by N 2 and O 2 
{p=Pe,N2 + Pe, 02 ) which Can be calculated by the following equations 

dPe,N2 


dr 


— ^N2{^) {Pat,N2 ~ Pe,N2] 


dpe,02 

dr 


— ^ 02 ( 1 ) {pat,02 - Pe,02) 


( 3 ) 


( 4 ) 


The photoionization rate of N 2 and O 2 is obtained by the Keldysh-PPT (Perelomov-Popov- 
Terent’ev) formula 361. 

The initial pulse investigated in this paper is Gaussian beam which can be written as: 


E{r, T, 0) = 


2P- 

' i/t 


ITTn 


p2 ^2 

exp ( —2 - 2 

^0 ^0 


( 5 ) 


where rp is the temporal duration. 
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III. RESULTS AND DISCUSSIONS 


Follow Li et al. 


17|, we define the collapse distance as the distance between the light 


source and the position where the laser beam has the smallest radius. Fig. [T] shows the 
evolution of the beam radius of pulses which have different temporal durations (80 fs and 
300 fs). It can be seen that the collapse distances of the two pulses have big difference (207 m 
and 135 m). In contrast, the semi-empirical formula in Eq. ([1]) gives = 129 m since 
it does not takes into account of the pulse’s duration. Therefore, the temporal duration of 
pulse may have a large influence on the collapse distance in some circumstances. 



FIG. 1. Beam radius of the pulse at different propagation distances. The input beam are Gaussian 
pulses with ro = 18 mm and Pin/Per = 20. 

From now on, we focus on the influence of the temporal duration to the collapse distance. 
In the simulations, we consider the Gaussian pulse with two different beam radius (ro = 2 mm 
and ro = 9 mm). 

Fig. [2] shows the variation of the collapse distances with different initial powers and 
temporal durations. Here ro=2 mm. Pin/Per varies from 1.5 to 6, and tq varies from 10 to 
100 fs. For comparisons, the collapse distances calculated by is also given and shown 

in solid line. We can observe that L*®™ can approximate Lc very well for the long durations 
such as ro=60fs, 80 fs and 100 fs, but it fails for the short ones, e.g., ro=10 fs, 15 fs and 
30 fs. In order to see more clearly how the collapse distances of the short-temporal-duration 
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pulses change with the input power, we show the results in Fig. [3] for ro=10 fs, 15 fs, 30 fs 
and Pin/Per being in the range between 1.5 and 12. It can be seen from this hgure that the 
collapse distances of pulses with these short durations have a maxima and tend to in 

the limit of Pin/Per going to inhnity. We dehne the input peak power which leads to the 
maximum collapse distance as Pmc- If Pin < Pmc, fhe collapse distance increases with Pin, 
and if Pin > Pmc, the collapse distance decreases with Pin. 



FIG. 2. Variation of the collapse distance of pulses with the initial peak power. The input beam 
are Gaussian pulses with rg = 2 mm. 


Fig. m shows the collapse distances for the pulses with ro= 9 mm. Pin/Per varies from 
2 to 30, and Tq =40 fs, 80 fs, 200 fs, 450 fs, respectively. From Fig. 0] we can get the 
same conclusion as shown in Fig. [2] and [3] that Lc has a maximal value when the temporal 
duration is very short. It is worth pointing out that “very short” is relative to the radius 
of pulses. For example, Tq= 80 fs is “very short” when ro=9 mm, whereas it is not “very 
short” when ro=2 mm. 

The phenomena that can not £t Lc well for the “very short” pulses is due to GVD. 
The normal GVD enables power exchanges between different time slices of the pulse and 


disperses the latter in time, and 
at smaller intensity levels 


37 


his may contribute to maintaining the pulse self-guiding 


39| . If the pulse’s temporal duration is long enough, the 


effects of GVD can be neglected when the propagation distance is less than the collapse 
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FIG. 3. Variation of the collapse distance of the short pulses with the initial peak power. The 
input beam are Gaussian pulses with rg = 2 mm. 


distance, and thus L*®™* can fit the collapse distance very well. However, if the pulse’s 
temporal duration is very short, the GVD effect will have a strong competition to the Kerr 
self-focusing effect when the initial power is small. Let’s consider an extremely short pulse 
with small initial peak power. In this extreme situation, the GVD effect overwhelms the 
Kerr self-focusing, and the collapse distance equals the propagation distance at which the 
peak intensity decreases to Per- At the same time, the temporal duration varies with z as 
Ti{z) = ro[l -|- {z/2LavDyY^'^ 40|, here Lqvd = 'ro/2/c" denotes the dispersion length 12|, 
and thus the peak power can be written as P{z) = Pin/Ti{z). The collapse distance in this 
case can be approximated by 


' extreme 


= 2L 


GVD 




( 6 ) 


This equation may account for the phenomena that the collapse distances of short pulses 
increase with Pin/Per when Pin < Pmc (Figs. EH), though Eq. |6]can not be used to estimate 
Le accurately in general cases. 

With closer inspections for Figs. EH we hnd that Le is proportional to Pin/Per when 
Pin < Pmc, and the slope is dependent on the temporal duration and the radius of pulse. 
Taking into account the GVD effects, we hnally obtain the new semi-empirical formulas as 
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FIG. 4. Variation of the collapse distance with the initial peak power {Pin/Per range from 2 to 
30). The input beam are Gaussian pulses with ro = 9 mm. The temporal durations of the pulses 
are 40 fs, 80 fs, 200 fs, 450 fs. 

follow 

^serm,ne. ^ i 0.225 Lgvd^^ + 0.320LnF , if P.n < Pmc , (7a) 

[ Lf™ X 1.878^ , otherwise , (7b) 

where Lj:,p = im^rQ/XQ is the Rayleigh length, which account for the contribution from the 
diffraction. is a dimensionless parameter and 

T semi 

N = , (8) 

^GVD 

which characterizes the relative importance of the Kerr self-focusing and the GVD. When 
N 1 GVD dominates, while for N <^1 Kerr self-focusing dominates. 

From Eq. (7), we can see that when Pj„ < Pmc the collapse distance is mainly determined 
by the GVD and the diffraction, otherwise, all the GVD, the diffraction and the Kerr self- 
focusing plays important roles. The input power threshold Pmc can be obtained via equaling 
Eq. (|7a)) and Eq. (iTbl) . The relations between Pmc, the duration tq and the radius tq of 
the laser pulse are shown in Fig. [5l It can be seen that Pmc increases with tq. At the same 
time, Pmc decreases with tq, and it reduces to Per in the the limit of Tq —?■ oo. 

In real applications of the new semi-empirical formula, we do not need the value of Pmc- 
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FIG. 5. Dependence of Pmc on the pulse’s duration and pulse’s radius. 


Notice that Eq. fl7a|) is a monotonically increasing fnnction of while Eq. fl7bj) is a mono- 
tonically decreasing fnnction of and the value of input power corresponding to the inter¬ 
section of these two functions is Pmc- Therefore, we can directly use min{Eq. (17a|) . Eq. (ItSI) } 
to obtain the collapse distance Lc for a given Pin, without knowing the value of Pmc- 

Fig. [6] shows the comparisons of the collapse distances between the new semi-empirical 
formula, the direct numerical simulations, and the conventional semi-empirical formula. It 
can be seen that the new semi-empirical formula agree with the simulation results very well, 
for both the long and short pulses. 

Now we turn to the case shown in Fig. 1, in which the collapse distances are obtained 
by numerical simulations. Basing on the new semi-empirical formula given in Eq. (7), we 
obtain the collapse distance of 214 m (207 m in simulation) for 80 fs and 134 m (135 m in 
simulation) for 300 fs. 

Eq. (7) can also be used to calculate the collapse distance at high pressures. Lgvd is 
inversely proportional to the pressure, and play more important roles at high pressure. For 
example, in the work of [l^, the authors obtained the collapse distance of the pulse with 
tq = 50 fs and tq = 1.2 mm, Pj„ = dP^j. is 2.3 m at 10 atm, while the new semi-empirical 
formula gives = 2.2 m, perfectly matched. 
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FIG. 6. The collapse distance predicted by the new semi-empirical formula. The input beam are 
Gaussian pulses with (a) ro=2 mm, (b) ro=9 mm. 

IV. CONCLUSION 

GVD plays important roles in the evolution of the laser pulses in air. For the short 
pulses, due to the competition between GVD and Kerr self-focusing, there exists a threshold 
Pmc for the initial power, with which the collapse distance has a maximum value. If the 
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initial power is less than the threshold, the collapse distance is proportional to the input 
power. If the initial power is larger than the threshold, the collapse distance decreases with 
the the input power, and is larger than that obtained by the conventional semi-empirical 
formula. Taking into account the effects of GVD, we present new semi-empirical formulas, 
which match the numerical simulations very well. The formulas can also be applicable to 
the cases for high pressures in which the GVD effects become more important. 
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